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We discuss gravitational perturbations in the Randall-Sundrum two branes model with radius 
stabilization. Following the idea by Goldberger and Wise for the radius stabilization, we introduce 
a scalar field which has potentials localized on the branes in addition to a bulk potential. In our 
previous paper we discussed gravitational perturbations induced by static, spherically symmetric 
and nonrelativistic matter distribution on the branes under the condition that the values of the 
scalar field on the respective branes cannot fluctuate due to its extremely narrow brane potentials. 
fT") ■ We call this case the strong coupling limit. Our concern in this paper is to generalize our previous 

analysis relaxing the limitation of taking the strong coupling limit. We find that new corrections in 
metric perturbations due to relaxing the strong coupling limit enhance the deviation from the 4D 
Einstein gravity only in some exceptional cases. In the case that matter fields reside on the negative 
' tension brane, the stabilized radion mass becomes very small when the new correction becomes 
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I. INTRODUCTION 



Understanding the role of extra dimensions has long been a focus of research. Recent developments in string theory 
stimulate a new possibility in a way to realize our Universe, i.e., the "braneworld," The braneworld scenario suggests 
that our Universe is realized on a brane embedded in a higher dimensional spacetime [HIS El (see also 0,0,0)- 

The explicit models introduced by Randall and Sundrum (RS) are simple but have attractive features 0, (jj ■ The 
two branes model, which was proposed earlier, is constructed by orbifold compactification of the 5D anti-de Sitter 
(AdS) spacetime, in which the two flat branes are on the Z2 symmetric fixed points. The single brane model, which 
was proposed later, is obtained by pushing the second brane to infinity in the first model. 
Oh, Since the bulk geometry of these models is warped, the behavior of gravi ty is not so trivial, and various interesting 
aspects of gravity in these models have been discussed [7I liil. [ToL ITTl ITa . IT3L 111 . fl5| . One fundamental but remarkable 
fact is that in the RS single brane model 4D linearized Einstein gravity on the brane is derived from 5D Einstein 
gravity 0,0,0- Also in the RS two branes model, 4D linearized Einstein gravity is restor ed [lfl 113 if the distance 
between branes, which we refer to as radius, is stabilized [iH ITlIll lllll ^ lMlHl^l^l^l2Dt l27tl28| . Although 4D 
Einstein gravity is approximately recovered in both models, the corrections arise in different manners. In the single 
brane model, the leading correction to Newtonian potential appears in the form of an inverse cubic potential, and 
therefore the correction is long ranged. This is because the mass spectrum of the Kaluza-Klein modes is continuous 
starting with (mass) 2 = 0. On the other hand, in the two branes model, the mass spectrum is discrete, and hence the 
correction becomes short ranged. As was pointed out in Ref. |l3|, the leading correction to the metric perturbation 
can be relatively large in amplitude but it is proportional to the local energy density of the distributed matter on the 
branes. 

The gravity beyond linear perturbations is also an interesting subject to study. For the single brane model, to 
study non-perturbative aspects of gr avity , many authors have discussed black holes in the braneworld |23, CM |3l|, l32l 
l33l l34l l35l l3rj l37l l38l 13^ . Hol l4lT lia . l43j. However, any black hole solutions that can be thought of as a state after 
gravitational collapse on the brane have not been found yet. On the other hand, there are studies of compact star on 
the brane |44J, ll5| . A pioneering work of numerically solving a relativistic star on the brane was done by Wiseman 
|46j . Another direction of research is to study higher order perturbations. Second order perturbations have been 



studied, and 4D Einstein gravity was proven to be restored under certain restrictions 
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For the two branes model, the recovery of 4D Einstein gravity is also concluded in our previous paper |49j under 
some assumptions. However, the mechanism for the recovery in each model is not so clearly understood as in the case 
of linear perturbations. 

The basic assumptions taken in the second order perturbations of the two branes model are that the radius is 
stabilized by the mechanism proposed by Goldberger and Wise (GW) [ly, and that the matter distribution on the 
brane is static, spherically symmetric and non-relativistic. In the GW model for radius stabilization, a bulk scalar 
field is introduced. This scalar field has potentials localized on the branes as well as a bulk potential. A further 
technical assumption taken in the previous analyses [lCt |49j is that the values of the scalar field on the respective 
branes are stuck to fixed values due to extremely narrow brane potentials. We refer to this simplified setup as the 
strong coupling limit. In this paper we generalize our previous analyses of linear and second order perturbations [49| 
relaxing the restriction of the strong coupling limit |51|. and study whether there arise observable effects and/or no 
pathological feature in the metric perturbations. 

The paper is organized as follows. In Sec. [H] we briefly review the formulation developed in Ref. |49j, which we call 
paper I, summarizing the notations and the basic equations. We do not give the explicit form of all the necessary 
formulas that are already presented in paper I to avoid repetitions of rather lengthy expressions. We quote the 
equations in paper I such as (1-1.1). Throughout this paper, we use the same notations as those in paper I except for 
the subscripts "pse" and "A," which are introduced in Eq. i|13|) . In Sec. IIHI we study the corrections caused by the 
fmitcness of the coupling strength. Section Hvl is devoted to the summary. 



II. BRIEF REVIEW OF THE FORMULATION 



In this section, we briefly review the formalism and the results presented in paper I. We consider second order 
perturbations in the RS two branes model assuming that matter distribution is confined on one of the branes and it 
is static and spherically symmetric. 

A. Notation and assumption 

The Lagrangian for the scalar field (p introduced for the radius stabilization is 

£ = -b ab v.av,b - v B &) - E v w&) s (y - v*) > (!) 

where Vb and V(±) are the bulk potential and the potential on each brane, respectively. 
In analyzing metric perturbations in the bulk, we use the "Newton gauge," in which 

ds 2 = e 2Y dy 2 + a 2 [-e A ^dt 2 + e B ^dr 2 + e c ^V W] , (2) 

where the metric functions A, B, C, Y, and ip depend only on y and r, and a(y) is the warp factor that is determined 
by solving the background equations (1-2.10). The metric functions and the scalar field are expanded up to the second 
order as 

A{r,y) = A^(r,y)+A^(r,y), (3) 
<p{r,y) = My) + ^ {1 Hr,y) + ^ {2) {r,y), (4) 

where (j>o represents the background scalar field configuration, which depends only on y. The metric functions A, 
B, and C are related to each other by the relations (1-2.5) so that they compose the transverse-traceless part at the 
linear level. 

Our interest is the gravity induced by non-relativistic matter fields confined on the positive and the negative tension 
branes, which are located at y = y+ and y = y_(> y+), respectively. The energy-momentum tensors of these fields 
are assumed to be given in the perfect fluid form as 

T±/ = a ± A diag{-p±, P±, P±, P±}. (5) 

The warp factor a± := a(y±) in the definition of the energy-momentum tensors is incorporated so that p and P 
become the physical quantities measured by using the metric on the positive tension brane. 
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The 5D Einstein equations give four independent perturbation equations. The two equations are the constraint 
equations for V (J) and </? (J) , which relate them to [Eqs. (1-2.12) and (1-2.13)]: 



-J— 5,( a 2 r (J) )^ 

2k0o<2 2 



Lp ( J) (r, y) 



2kg 



d y A 1 S,p 



(6) 
(7) 



where we have introduced defined by e^ 1 ^ = and e*- 2 -* = 1. Hence, as for the scalar type perturbation, once we 
solve the perturbation Y^ J \ other variables and ip^ are also known. The other two equations are the master 
equations for and Y (J *> [Eqs. (1-2.14) and (1-2.15)]: 



a 2 a 2 



(a 2 A^) 



Sa ■ 



(a 2 Y^) = e^S, 



(8) 
(9) 



To solve the master equations, we must specify the boundary conditions on the branes. The boundary condition 
for A^ is given by Israel's junction condition, whereas the boundary condition for the scalar type perturbation is 
derived by integrating the equation of motion for the scalar field across the branes. It is well known that these 
junction conditions are easily obtained in Gaussian normal coordinates x a in which a brane is located at y — const 
hypersurface. Here we associate an over-bar with quantities written in Gaussian normal coordinates. The boundary 
conditions in the Newton gauge are derived by applying infinitesimal gauge transformations x a = x a + £ a to those 
written in Gaussian normal coordinates. These transformations are described in Sec. II B in paper I. 

The boundary condition for the scalar type perturbation is given by [Eq. (1-2.37)] 



'■P 



(-/) 



U) X 



2 1 jun 



(at y = y±), 



where we have introduced a coupling constant 

A± := 



(10) 



(11) 



An over-dot denotes differentiation with respect to y, and (f>o± := 4>o(y±). Note that the explicit expression for Sj un 
is given in Eq. (1-2.39). The parameter of the gauge transformation in y direction, £±( r )> is a function of r, and it is 
related to the energy-momentum tensor on the corresponding brane by Eq. (1-2.34). 

In paper I, the problem was simplified by taking the strong coupling limit. The condition of the strong coupling is 
that V^_) is sufficiently large. Taking this limit, we obtain 



Ah 



0. 



(12) 



In this paper we discuss the effects due to the terms that arise by considering non- vanishing A±. For brevity, we refer 
to these terms as interaction terms. 



B. Derivative expansion 



Formal solutions for perturbation equations JBJ) and I® were derived in paper I. The transverse-traceless per- 
turbations A are decomposed to massless zero mode Aq and massive mode As- Using the zero mode truncation 
approximation, the long-ranged part of the transverse-traceless perturbation Aq is evaluated [Eq. (1-3.5)]. The re- 
maining part A$ [Eq. (1-3.9) and (1-3.12)], which arises due to the contribution from massive Kaluza-Klein modes, 
is evaluated by using a derivative expansion method. In this method, we expand perturbation variables in terms of 
the expansion parameter (Hr^)^ 1 ~ H~ 1 d r assuming that the typical length scale r+ of perturbations is much longer 
than the 5D curvature scale H~ l := (d/a) -1 . It is important to stress that this derivative expansion method is valid 
only when the mass of the first excited mode is sufficiently large. In the limit ?/_ — ► oo, the excited mass spectrum 
becomes continuous, and therefore the derivative expansion method is no longer valid. (See Sees. Ill B.l and IV B in 
paper I for more details.) 
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As for the scalar type perturbation, there is no zero mode owing to the stabilization mechanism ^(| • To discuss the 
contributions from massive modes, we expand the perturbation variables by using the derivative expansion. Although 
the massive modes seem to give only the short-ranged part, it turns out that the formal solution of the lowest order 
in the expansion includes long-ranged metric perturbations. Hence, we refer to this part as the pseudo long-ranged 
part. The formal solution of the next order, which is obtained by an iteration, is referred to as the short-ranged part 
Ys- (See Sec. Ill C in paper I.) 

The pseudo-long-ranged part includes the contribution from the interaction terms that are higher order in deriva- 
tive expansion as the short-ranged part. Hence we further divide the pseudo-long-ranged part into two pieces; the 
contribution from the interaction terms, which we denote Y\, and the remaining terms, which we denote Y pse . We 
quote the explicit expression for the pseudo-long-ranged part [see Eqs. (1-3.21) and (I-B.6)], 



y^+Yl J) +Yy\ 



AF p (/ e )(r,j/±) = - 



Ay A (J) (r,y±) 




dy(a 2 v±A\ 



V H ^ASy — a 2 S^ 

2n<j) 2 



(13) 



where 



Li J Xr) = H{y a ) 



e^a 2 f 



a \ oi. 3 un 



<?W 2 3<^AF \ 



20g 



2ko 2 0q 



Here N is the normalization factor defined by 



N~ 



:= 2 



a dy, 



v+ 



and u± is given by 



u±{y) := 1 - 2Hv ± , v±(y) := — / a 2 dy' 



y f 



(14) 



(15) 



(16) 



We mention that the source term Sj un that is defined in Eq. (|1(J|> contains V/±\ an d V/±y and hence X±Sj : un does not 
vanish even in the limit (|12f> . However, it was shown that the combination that appears in Eq. (|14|) vanishes in the 
strong coupling limit. Namely, we have 



L 



(J) 



(for A± 0) . 



(17) 



Note also that, in Eq. (|13|l . the number of r derivatives on the right hand side is larger by two than that on the 
left hand side. This fact is manifest for linear perturbations, and it also turns out to be the case for the second 
order. Hence, when we evaluate Y^ iteratively, the leading term with respect to gradient expansion is obtained by 
substituting pseudo-long-ranged part Y^e into Y^ on the right hand side of Eq. (|T4*jl . Although the short-ranged 

part Yg 1 ^ also contains interaction terms, we do not discuss them in this paper because they are even higher order in 
rj. 

At the linear level, the transverse traceless part composed of A, B, and C does not have interaction terms. The 

(i) 

parameter of the gauge transformation is also unaltered by the effect of non- vanishing X±. On the other hand, 
(pt 1 } is related to and changes according to the change in Y"W . 

By repeating the derivation of the expression for the temporal component of the metric perturbation induced on the 
branes (1-5.1) starting with the equation for gauge transformation (1-2.25), we will find Eq. (1-5.1) is slightly modified 
as 



AA { J \r,y ± ) = 8irG(p± + 3P±') - A i v ± A±^ y + £±Ab,r + #(£±) 2 +2JV£V<4(S£ - S£) 



v- 

2 V / dy 

'y+ 



o, 2 (Sa± - S^±) + a 2 v±AS v ± 



3uj 



-AS 



Y± 



AY, 



(J) 



(18) 
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For the spatial components, it is convenient to take the isotropic gauge. We simply quote Eq. (1-3.15) in which the 
isotropic gauge is taken in Jth order: 

AB«(r,y ± ) = - l -AA^ T - \aY^ | e O + + 2Ha%Sf) . (19) 

The explicit additional term in these equations is only the last term in Eq. (|18fl . but there are implicit changes through 
Y~( J ), (p^ J \ and A^ J \ The explicit expressions for the source terms S* and §* are given in Eqs. (1-2.16), (1-2.18), 
(1-2.19), (1-2.21), (1-2.33), (1-2.35) and (1-3.14). 



C. Corrections 

As we have done in paper I, we assume that matter fields reside on one of the two branes. By this simplification, 
the sum of the Newton potentials <I>±, which are defined by A$±(r) = nNp± /2, is replaced as 

5>*->*± (y = y±). (20) 

Then, the long-ranged part of the transverse traceless metric perturbations is given by (1-4.7) and (1-4.8) as 

[$&{r,V±)]sc = -[Bi l 2(r >y± )]sc = 24>± , (21) 

where [• • - ]sc means the quantity in the A± — > limit. The index ± associated with the metric functions (not with 
A and L) specifies the side on which matter fields are distributed. 

To obtain an approximate estimation for the short-ranged part, we assume that the back reaction of the bulk scalar 
field to the background geometry is weak; \H\/H 2 < 1. In this case the metric approximately takes the AdS form 

atone-to", (22) 

where £ is the curvature radius of the bulk. For later convenience, we evaluate iV in this weak back reaction case. It 
is approximately given by 

taking into account the fact that the integral l|15[) is dominated around y = y + . The second term in the round brackets 
is the leading order correction due to the back reaction. Hereafter we set a+ = 1. 
Under the assumption of weak back reaction, it was shown that [ToL l49| 

[4l(r,y ± )] SG =p ± 0(rlA$ ± ), (24) 

where the suppression factors are 

,_ = JL = (^)Vi^)V') 2 . (25) 

The short-ranged part for the scalar type perturbation is similarly suppressed as 

[Y^{r,V±)] sc = ^0(rlA^ ± ), (26) 

where m| is the mass squared of the so-called radion, which is the mode corresponding to the radius fluctuation, in the 
strong coupling limit. To be precise, the radion is defined by the mode with the lowest mass eigenvalue in the scalar 
type perturbation. Since m| < O(o 2 _£~ 2 ), which is given by Eq. I)29|) below with A± = 0, the corrections from the 
scalar type perturbation tend to be larger than those from the Kaluza-Klein modes in the tensor type perturbation. 

On the positive tension brane, the short-ranged part is suppressed when the typical length scale r* of perturbations 
is much larger than a^rrig 1 . On the negative tension brane, the short-ranged part is suppressed for > 0.1 (a_/mg£) 
mm, where the ratio a + /a_ is set to 10 16 , the value to solve the hierarchy between Planck and electroweak scales. One 
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may think that this effect is potentially observable especially on the negative tension brane. As we have mentioned 
in the introduction, however, the short-ranged part is proportional to the local matter energy density p± oc A$±. 
Therefore, the short-ranged part dose not contribute to the force outside the matter distribution, and it is not observed 
as a change of the Newton's law. Not the force but the change of the metric perturbation due to the short-ranged 
part becomes significant on the negative tension brane only when > (m|£ 2 /a?_)0(TeV 4 ). 

Extending the linear analysis mentioned above, second order perturbations were investigated in paper I. In the 
case that matter fields are confined on the negative tension brane, the long-ranged part at the second order contains 
corrections to 4D Einstein gravity with relative amplitude of 0(/3_ / a_) compared to the ordinary post-Newtonian 
terms. The relative amplitude of the corrections at the second order looks less suppressed by an extra factor of l/a 2 _ 
than that at the linear order. However, in the strong coupling limit, it was shown that these enhanced corrections 
are completely canceled by the contributions from the short-ranged part, and the recovery of 4D Einstein gravity was 
confirmed. In the case that matter fields are confined on the positive tension brane, the corrections in second order 
perturbations are relatively 0(/3+/a?_) compared to the usual post-Newtonian terms. However, the appearance of 
the enhancement by the factor of l/a?_ is very likely to be an artifact due to the gradient expansion method. Since 
the condition that the typical length scale of spatial gradient is larger than that of the change in the fifth direction 
becomes (£ 2 /a 2 _r 2 ) = (/?+/a 2 „) <C 1 near the negative tension brane, (3 + /a 2 _ appears as an expansion parameter. 



III. CONTRIBUTIONS DUE TO NON- VANISHING A± 



In the analysis taking the strong coupling limit, we neglected the terms that vanish when A± is set to zero. Here 
we consider the effect of non- vanishing X± to remove this technical limitation. We first discuss linear perturbations, 
and after that we study second order perturbations. 



A. Mass spectrum 

Before discussing metric perturbations induced by matter fields, we study the change of mass eigenvalues for the 
scalar type perturbations due to the effect of A±. Setting £± = in Eq. (1-2.41) with J = 1, we obtain an eigenvalue 
equation for Y"W as 



1 



2 k 



a 2 <t>ldy—^d y a 2 - —a 2 4>l + m 2 



(\ + 2^A CT( %-y ff )^ F« 



0, 



(27) 



where the 4D Laplacian operator A was replaced with the mass squared m 2 . A general solution for small m 2 ignoring 
boundary conditions at y — y± was approximately constructed in Ref. .IQj as 



1 _9 ~ ' 



3(j> 



a 2 cf> 2 



j 2k m 2 



dy' + c 



(28) 



where c is a constant and J\f is a normalization constant. Imposing the boundary condition I|1U|) on this approximate 
solution, we find that c = — X+m 2 / (f>Q + a^_ and a mass eigenvalue close to zero is given by 



2k 

T 



r,2 



dy' 

a 4 j> 2 ' ^ <■'■'■- 



(29) 



Using Eq. (I-B22) with the assumption of weak back reaction, the above expression is approximately rewritten as 

-l 



™2 ^ 2 



4 / 2 



a 2 _N 2 



where we have introduced non-dimensional parameters related to the coupling A± by 

27V 3 A± 2A^ 3 A± 



a± := 



k^I ± 



3H± 



(30) 



(31) 
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When A + > nij 2 ! -1 , the expansion (|28() is no longer valid since the correction that comes from the constant c 
becomes larger than J\f/a 2 . In such cases, instead, we can consider the large A + limit keeping m 2 A + finite. Since we 
can neglect to 2 term in the bulk when to 2 <C trig, we find that an approximate solution to the above equation is given 



by Y"W w U+{y) hi this limit. From the condition l|10|) at the boundary y 
this mode is determined as 



y + the mass eigenvalue corresponding to 



m 2 = 



2H+v + (y + ) 



2 

a7? 



(32) 



where we used Eq. (I-B3) with the aid of Eq. IL'-it . 

When A+ (or A_) takes a large negative value, we can see that the above mode of small mass becomes tachyonic, 
and hence such a model is manifestly unstable. As we decrease |A±| starting with A± = — oo, the absolute value of m 2 
in the expression (|3L)|) increases from m 2 = and diverges to | — oo| when A± is a certain negative value A±, which 
depends on the details of stabilization model, e.g., <j>o±- Then the mass eigenvalue returns from +oo to m| as X± 
increases from A±. However, this does not directly indicate that the tachyonic mode disappears for A± > A±. This 



is because Eq. lKi()[l is no longer valid for relatively large \m 2 \ 
the WKB approximation as 



For large v 



-m 2 , one can solve Eq. (|2"7Jl by using 



y(l) Q ie " 'f V dy'/a{y') _|_ (J 2e -v f y dy' /a(y') 



(33) 



where C\ and Ci are constants. Imposing the boundary conditions (|1UI) at y = y±, we obtain (X+v — a+)(A_f — 
a_) — exp[2z/ a _1 dy](A+^ + a + )(A_^ + a_) « 0. Without loss of generality we can assume that v is positive. Then 

the exponential factor exp[2z/ a~ 1 dy] is very large, and hence a solution to the above equation is approximately 



obtained when either A+ or A_ is negative as 



max(— a+A + 1 , — a_A_ 1 ). 



(34) 



This eigenmode is the anticipated tachyonic mode, which remains to exist for any small negative value of A+ or A_. 
Although the values of A± depend on the details of the stabilization model, it is natural to consider the case in which 
A±| is less than or equal to 0(£), and A± must be positive. To conclude, we find that the model has tachyonic mode 
if A+ or A_ is negative. 



B. Linear perturbation 
From Eq. Ijl8|l . the temporal component of the linear perturbation is given by 



(35) 



As for the spatial part, we obtain from Eqs. (|19|l and (|35|) 



(36) 



Substituting Eq. Q into Eq. O, we can evaluate L± by iteration. Keeping up to the order r+ 2 , we obtain 



a± 



(37) 



Then Eq. lfT^)l is evaluated by using Eq. as 



Y^(r,y-) w 2 A 



Y$(r,y+) « 2 A 



' H+ + H_ 



?>H%J H 4 



a\H_ 



= /3 + 0(r 2 A$ + ) 



= /3_0(r 2 A$_) [a 4 _a+ + a-] 
H+ \ 2 



3Hj 



(38) 
(39) 



On the negative tension brane, this new correction becomes important compared to that from the short-ranged part 
when the factor in the square brackets on the right hand side exceeds 0(a 2 _/m 2 s £ 2 ). In particular, when a^a + 3> 1, 
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Eq. applies and we find m 2 w 47V 4 /(3a+|ff+|) < atH\/(l 2 \H + \). Because of the factor at, the mass of the 

stabilized radion becomes even smaller. On the positive tension brane, the first term in the square brackets on the 
right hand side is suppressed only by the factor H + /H^_, which is small but is not hierarchically suppressed. As 
long as A+ takes the natural order of magnitude smaller than £, a+ is at most 0(H+/H+). Then, the correction to 
Y\+{r, y + ) stays less than 0(j3 + )0(r 2 A$ + ). However, when A + is much larger than £, the correction becomes larger 
by the factor of X+/£ than that in the strong coupling limit. Although these choices of parameters are not natural, 
the possibility of the enhanced correction without changing the order of radion mass might be interesting. 



C. Second order perturbation 



1. Temporal component 



Let us discuss the contribution due to interaction terms to second order perturbations. In the strong coupling limit, 
the leading terms in second order perturbations are shown to be identical to that given by 4D Einstein gravity; 

[A4 2 2] sc = ^G{pf + 3Pf ) - 4$ ± A$ ± + O [y^) ■ (40) 

We will show that the corrections due to interaction terms are similarly suppressed as 0(\/a*_a±rt), where and 
hereafter we assume that X± is not hierarchically enhanced and hence a± is at most 0(1). In the following discussion 
we concentrate on the terms of 0(r~ 4 ) with respect to the derivative expansion, neglecting the higher order terms 

(i) 

than 0(r~ 6 ). For simplicity, we adopt as a choice of radial gauge in linear perturbations, keeping second order 
perturbations still in the isotropic gauge. 

We quote the dependence of the first order perturbation variables on the warp factor from Eq. (1-5.4): 

A {1) ~n° Y {1) w {1) + 1 P> ~ — f4Fl 



As for u± and v±, we have Eqs. (1-4.11) and (1-4.12) as 



2 2 
at 



u± ~ -f + 1 , v ± ~ -f - 1 . (42) 



or a 



We list the perturbation variables of the first order that have a correction due to non- vanishing A± . As we have 
discussed, there is correction to which is denoted by in Eq. i|13|) . Since ip^ is related to by Eq. (1-2.13), 
ip^ also has correction accordingly, which we denote by ip^ . In the following discussion, the values of tp^ evaluated 
at y = y± are necessary. They are obtained from Eq. iJTUj l as 



2na±<p ± V r */ 



(43) 



Although the source terms for second order perturbations are mostly written in terms of the variables in the Newton 
gauge, the expression A,y in Gaussian normal coordinates is also necessary. From Eq. (1-2.31), we can read the 
interaction term in A$ as 

d y A?=T 2 -^A 1] (aty = y± ), (44) 



where we used tp\\y±) = tp\(y±), which follows from the fact that the gauge transformation (1-2.26) of tpW is not 

altered by the effect of non- vanishing 
From Eqs. (|13|l and 143(1 . we obtain 



altered by the effect of non- vanishing A±. Also for the variables Y$ and <p^\ we give the dependence on a_ and r + 



Y^(y ± ), ^(tfi)--^. (45) 
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For later use, we quote the relations (I-B12) and (I-B13) in the notation of present paper as 



y« + 2H ± I 



(1) 
y 
± 



J2 (at y = y±) 



<T=± 



Y. 



(i) 



2^± (i) 



®* (at y = y±) 



(46) 



<T=± 



We begin with the case that matter fields are confined on the negative tension brane. We discuss the contributions 
from each term in Eq. ljl 8|) one by one. The second term on the right hand side of Eq. Ijl8|l gives 



y \2 



(47) 



The contribution from the source terms and Ss is evaluated by using the expressions given in Eqs. (1-2.33) and 
(1-2.35) as 



27V ™t(SI - S£) = -2Nat [5f -S s ] sc + 0( —^-Y\- (y 



O 



(48) 



As for (Sa — S,[,), the underlined terms in Eqs. (1-2.16) and (1-2.19) give the leading correction as 



2N 



a 2 {S A 



' ■ \sc 



O 



y+ 



^y A _( y _) ) o 



(49) 



where we have used the fact that <p\-(y-)/<f)Q- — 0(.ff~ Yx_(y_)). The leading order correction from S v comes from 
the last two terms in Eq. (1-2.18). These terms are rewritten as in Eq. (1-5.7), in which we did not assume the strong 
coupling limit. The underlined terms on the right hand side of Eq. (1-5.7) give the leading order correction as 



2N / dy a'v-ASy 
'v+ 



2 k 



JSC 



- 2NA / dy 



a 2 u. 



(Y 2 -Y 2 ) 

V 1 psc/ 



+o 



Yx-(y-] 



o 



The last term cancels the contribution from Sy that is given by 

AS Y = [AS Y ]sc + ^a 2 <p 2 A(Y 2 - Y 2 SC ) + O 



(50) 



(51) 



where again the leading correction comes from the underlined terms in Eq. (1-2.21). The second term in Eq. I|47|l and 
that in Eq. I|5U|) . which potentially give enhanced correction to of 0(£ 2 <J? 2 _/r 2 a'l), cancel each other with the aid 
of Eq. (1-3.27). The terms of O (r~ 2 $_Y\- (?/-)) are smaller by the factor of <3?_ than the correction we have found 
for linear perturbations. 



Now we consider the last term of Eq. I|18fl . i-e., the contribution from Yc ' . To evaluate F A , we study L_ given 
in Eq. Q14|). According to the dependence l|41|l . some terms have possibility to give an enhanced correction to AY$ 2 ' 

of order 0(l/a^_rf). Note that the terms of 0(1 /atxt) m AL_ give the terms of 0(l/a1rf) in AY^ 2 \ while the 
contribution from AL + does not change its order with respect to a_. In the following discussion, we keep only the 
relevant terms that might give enhanced contributions of 0(l/a1r* ) to AY {2) . 
Keeping the terms of 0(1/ 'a^r 2 ), the source term Sj un in Eq. fl4^l becomes 



A- 4>Qip 2 f 3ipAY 

2o„ jun 2o;j 



2Ka 2 _<j>Q 



(e_,) 2 

2a 2 



2Ka 2 _4>Q_ 



Y + a 2 _HY' 



(52) 



J y=y- 



where we used the expression given in Eq. (I-B18). Here the term of §y is evaluated by using Eq. (1-2.21). Keeping 
the terms of 0(l/alr 2 ), we obtain 



§ Y + a_HY 2 



3Y r AY + ^-^{AY)^ r + —((p. r A(p - <p(A<p) tr ) + 4a 2 B^ (2HY r + H^- 
0o 3 ' V O 



dr , (53) 
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where we have used Eq. (|46|l . The right hand side of Eq. I|14fl also contains Y^ 2 \ To obtain the lowest order 
correction, this Y^ can be replaced with Ypsc ■ We can neglect the contribution from Yp S 2 2 (r, y+ ) in Eq. I|13|) . which 
is not enhanced with respect to the factor of a_. The relevant terms in AFp 2 i(r, y„) are O^l/at^r 2 ), which are 



Ai'p'2 (r,y. 



S^j + 2Ha 2 _S^ 



(54) 



where the first term comes from the integration of S v . Using Eq. (1-2.16), S^, in the second term can be explicitly 
written down. The relevant contribution comes from only the underlined terms in Eq. (1-2.16). With the aid of 
Eqs. (1-2.10), (1-2.13), (I-B14) and the fact that (a 4 5 a ) ia = 0(a 2 /r 2 ), we obtain 



dr\ Aa 2 B ty (HY r + H 



As for , we can read from Eq. (1-2.35) as 



<P.r 
00 



r8/3 



9, 



,8/3 



( 3 -Y 2 

\2 ' r 



a 2 _s; 



a±B iV Y r 



e(AF), r dr-eAY-H(e__, r )\ 



where we have again used Eqs. (|46|l . Substituting all the results into Eq. I|14fl . we obtain 



3a- 
4/V3 



dr 



K i 2 

3 ( ^ 



• ~ 2 2ff(AF), . ~ 

»05,r) ),rH • (^-<PoC 



r 



(55) 



(56) 



(57) 



Further application of Eqs. I|43|l and (|46|l reduces the order with respect to either a_ x or r*. Therefore we find 

LL 2) =a-0(^ 2 $ 2 _/r2 a 2 _) + 0(l/r4). Although we have not discussed the contributions associated with the factor a+ 
in detail, it is manifest that they do not have any enhancement with respect to the hierarchy factor of l/a_. Therefore 
the interaction terms in second order perturbations become 



A4 2 J(r,2/_)= \AA™(r,y-) 



O 



sc 



[0{ata+) + 0(a-)] 



(58) 



which are suppressed compared to the correction at the linear order by the factor of 

Next we consider the case that the matter distribution is concentrated on the positive tension brane. In this case, 
the first order quantities listed in Eqs. (|41|) and Q42JI do not suffer enhancement with respect to the factor of l/a_. 
The factor l/a_ arises only through Y\+(y-) and (fx+(y-), both of which are 0(l/a 2 _). Note also that the terms 
with A + are always associated with u + (y + ) = 0(H + /HV). From these observations, we find 



O 



sc 



(59) 



The amplitude of this second order correction is not simply suppressed by the factor of compared to the correction 
at the linear order. There is a difference with respect to the power of a_ in the term with In paper I, we have met 
a similar phenomenon in the analysis of correction due to Kaluza-Klein modes at the second order. As was discussed 
there, this phenomenon is a natural consequence of our gradient expansion approximation. 



2. Spatial component 

Finally, we comment on the spatial component. The last two terms on the right hand side of Eq. i|19|) are evaluated 
in the same manner as was done for AAq ^ . The source terms except for Sb have been already computed in evaluating 
A A a ± . The contribution from Sb is given by 

AS B w [AS B ] sc ~ 4oA[el(^ - Ml)] , (60) 
which is similar to that given in Eq. (|47|l . Combining all, we obtain 

B=[B] SC + 0(A ( 2 1-[A \% C ) , (61) 

and therefore the correction is the same order as that of the temporal component. 
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IV. SUMMARY 

We have discussed metric perturbations in the Randall-Sundrum two branes model with radius stabilization. As 
a mechanism for radius stabilization, we have assumed a scalar field with a potential in the bulk and that on each 
brane. In our previous work (paper I), we took the strong coupling limit, in which the brane potential is extremely 
narrow so that the values of the scalar field on the branes cannot fluctuate. In this paper, we extended the previous 
analysis relaxing the limitation of taking the strong coupling limit. 

In the strong coupling limit, it is known that the mass squared of the stabilized radion tends to be hierarchically 
small as m 2 s = 0((a-/a + ) 2 £~ 2 ), where I is the bulk curvature scale and (a_/a+) is the ratio of the values of the warp 
factor on the respective branes. First we examined the shift of this mass eigenvalue when we relax the limitation of 
taking the strong coupling limit. We have shown that a tachyonic mode appears when either of the coupling constants 
A± defined in Eq. (|11J) is negative. Hence, the models with such parameters are unstable. When both A+ and A_ 
are positive, we derived formulas for the mass squared of the stabilized radion in Eq. (|32[l for A+ nig 2 ^ 1 and in 
Eq. (EOJ for A+ <C m^r 1 , The mass squared of the stabilized radion is affected by A_ when A_ > I while A+ only 
when A + > m s 2 l _1 . These are rather exceptional cases since the order of A± is typically less than or equal to O(i). 

Next, we have examined the effects on metric perturbations induced on the branes by matter fields up to second 
order, assuming that the matter distribution confined on the branes is static and spherically symmetric. For simplicity, 
we assumed that matter fields reside on either of the two branes. The results for the case that the matter fields are 
on the negative tension brane are summarized by Eqs. (|38|l and (|58() . where a± and (3± are defined in Eq. 131|) and 
Eq. (|25|l . respectively. The correction due to the finite coupling becomes important compared to that already existing 
in the strong coupling limit only when the mass of the stabilized radion is significantly reduced by the effect of 
nonvanishing A± . The results for the case that the matter fields reside on the positive tension brane are summarized 
by Eqs. (|3*§|) and (JSSJ. It is possible that the correction due to the finite coupling becomes important when A+ becomes 
much larger than £ without changing the mass of the stabilized radion. The corrections are enhanced by the factor of 
compared to those present in the strong coupling limit. The result for second order perturbations l|59ll seems to 
show that the correction associated with A_ are enhanced by a factor of a^/a/L. However, this is an artifact due to 
the limitation of the present approximation using gradient expansion method. 

In conclusion, under the assumption that the matter distribution is static and spherically symmetric, we have 
confirmed that 4D Einstein gravity is approximately recovered up to second order perturbations relaxing the limitation 
of taking the strong coupling limit. The condition that the corrections due to a±, (3±, and radion mass are sufficiently 
suppressed gives a consistency check for any stabilization model to use the scalar field of Eq. Q [5(| • 
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